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INTRODUCTION 


2 

A  new  method  of  finding  prime-rich  equations  of  the  type  I=x  -x+c  (c=  2N- 

1,  N=1,2,3 . )  has  been  devised.  The  Euler  type  equation  with  the  highest 

known  density  of  primes  is  I=x  -x+41  where  the  density  Dc  is  defined  as 
P(xN)100/xN  percent;  P(x^)  is  the  number  of  primes  and  xN  is  the  limiting 
value  for  the  range  x=l  to  x^.  Unless  otherwise  stated  Dc  will  be  expressed 
as  the  percentage  for  xN=40;  for  this  is  100%.  By  using  this  method 
fifteen  equations  (equivalently  fifteen  values  of  c)  were  found  which  have 
higher  densities  of  primes  than  when  xN  is  set  to  equal  to  10,000.  Table 
1,  shows  these  equations  plus  five  additional  ones  whose  densities  are  nearly 
equal  to  D^. 

PREVIOUS  WORK 

The  method  is  based  on  an  analysis  reported  previously.*  There  a  special 

diagonal  array  of  the  odd  integers  was  shown  to  be  correlated  with  equations 
2 

if  the  form  I-x  -*+c,  c  =1  for  the  first  column  of  the  array,  c=3  for  the 

2 

ond  etc.  Similarly  the  rows  were  described  by  x  +x-r  where  r=l,3,5,.... 
fo r  the  rows  N=l,?,3,...  respectively.  An  analysis  of  the  rows  is  a  separate 
■,*  id>  ’'i  i'self  hut  was  not  pursued  in  reference  1  and  is  not  considered 
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TABLE  1.  EQUATIONS  x2-x+c  FOR  c=41, 21377, .  WHERE  THE  DENSITY,  D  IN 

PERCENT  (100xPRIMES/40)  FOR  x=l  TO  40  IS  >  67.5%  AND  WHERE  Dc 


IS  ALSO 

GIVEN 

FOR  xN 

=100,  ; 

200,  ETC. 

Equations 

havi ng 

hi gher 

D  values 

then  c=41 

for  x^=  100000 

c 

xN=4° 

100 

200 

400 

1000 

2398 

5000  8000 

41 

100 

86 

78 

67.5 

58.1 

50 

45.22  42.575 

21377 

80 

75 

72 

66.75 

57.5 

51.33.. 

45.78  43.5 

204791 

77.5 

57 

58.5 

56.75 

54.7 

49.66.. 

45.68  43.0125 

27941 

75 

77 

74 

68.75 

59.9 

53.25.. 

48.06  45.625 

55661 

75 

73 

72 

69 

62.1 

54.92.. 

49.96  46.7875 

22697 

75 

69 

69.5 

65 

56 

49.66.. 

45.62  43.0125 

72491 

75 

68 

69 

67 

61 

54.37.. 

49.06  46.2875 

52727 

75 

66 

64 

60 

54.5 

49.87.. 

45.32.  42.5625 

74147 

75 

65 

64  5 

61.75 

56.3 

50.91.. 

46.38  43.8125 

41537 

70 

73 

71 

68 

60.6 

53.33. . 

48.78  45.8 

595937 

70 

60 

62.5 

63 

60.9 

57.96.. 

53.88  51.05 

47891 

67.5 

69 

67 

63.75 

57.2 

50.66. . 

46.02  43.4875 

878357 

67.5 

65 

60.5 

57.25 

53.6 

50.625. 

48.28.  46.1625 

579431 

67.5 

64 

57 

57.5 

55.7 

53.25.. 

49.18.  46.9375 

310691 

67.5 

62 

59.5 

60 

58 

53.66. . 

49.52.  47.225 

Equations  with  Dc  values  close  to  those  for  c=41 

19421 

75 

74 

70 

64 

55.8 

48.37.. 

44.12.  41.3625 

31511 

67.5 

66 

64 

58.75 

54 

48.49.. 

43.76.  41.2 

54407 

67.5 

66 

62 

61.25 

55.7 

50.25.. 

44.82.  43.325 

41201 

67.5 

65 

63.5 

60.5 

54.7 

48.74.. 

44.04.  41.55 

221717 

67.5 

57 

56.5 

54.5 

51.5 

47.58.. 

44.38.  42.25 

%  .*• 

vv>y^ 


10000 

41.49 

42.55 
41.82 
44.66 
45.44 
42.07 

45.33 

41.55 
42.85 
44.80 
49.78 

42.33 
45.10 
45.89 
45.91 

40.24 

40.32 

41.13 

40.24 

41.16 
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here.  From  the  special  diagonal  array  a  set  of  arrays  was  derived  called 
primitive  cell  arrays,  one  for  each  prime  value  of  c.  Each  such  array  has  the 
value  c  substituted  for  every  integer  in  the  original  diagonal  array  which  is 
=  0(mod  c);  at  each  integer,  I, not  =  0(mod  c)  the  space  (corresponding  to  its 
position  in  the  original  array)  is  left  blank.  Reference  1  included  tables 
showing  primitive  cell  arrays  for  3,7,11  and  a  modified  one  for  5. 

Other  arrays,  not  included  in  the  report,  were  constructed  for  c=P  (P  is 

prime)  =  13,17,19,23, . ,41.  It  was  noticed  that  in  each  primitive  cell 

array  there  were  "empty"  columns  (and  rows)  -  by  empty  is  meant  columns  in 
which  integers  were  =  0(mod  P).  The  number  of  empty  columns  can  be 
expressed  as  Nec=(P-l)/2  e.g.  for  P=3  Nec=l ;  for  P= 1 3  N0  cq1  =6.  It  was 
proved  that  for  the  primitive  cell  arrays  of  3  and  and  5  the  patterns 
are  repeated  throughout  these  arrays  no  matter  how  far  the  columns  and  rows 
are  extended  (for  any  P  the  primitive  cell  pattern  consists  of  a  PxP  "area" 

O 

made  up  P  spaces  occupied  by  P’s  and  P  -P  blank  spaces).  It  was  assumed,  with 
some  supporting  evidence,  that  this  holds  true  for  every  such  array.  Table  2 
summarizes  the  equations  which  describe  the  number  of  sets  of  empty  columns 

for  c  =  P  =  3 ,5 ,7 , 1 1 . . ,37.  For  the  PCA  of  3  (PCA  will  be  used  instead  of 

primitive  cell  array  for  the  sake  of  brevity)  the  empty  columns  consist  of 

just  one  set  Ngc=3N+3,  N=0,l,2,3 .  It  can  be  seen  that,  if  one  ignores 

the  presence  of  all  other  divisors,  a  column  of  this  set  can  have  a 
density  of  100?  e.g.  c=41.  In  contrast  those  columns  which  do  contain 
integers  which  are  congruent  to  0(mod  3)  are  described  by  the  two  sets  3N+1 
and  3N+2,  N=0,l,2,3 .  The  former  set  describes  the  1st,  4^,  7*h, . 


columns  of  the  cell  of  3  for  which  every  third  integer  of  every  column  is 


TABLE  2.  EQUATIONS  SHOWING  THE  NUMBER  OF  SETS  OF  EMPTY  COLUMNS  FOR 
EACH  PRIME  VALUE  OF  P  FOR  THE  PRIMITIVE  CELL  ARRAYS 


p 

No.  of 

Sets  of  empty  columns 

sets 

per  array  (N=0,1 ,2, . . ) 

3 

1 

3N+3  * 

5 

2 

5N+1 ,  5N+4 

7 

3 

7N+2,  7N+6,  7N+7 

11 

5 

11N+1 ,  11N+4,  11N+8,  11N+9,  11N+10 

13 

6 

13N+2,  13N+3,  11N+8,  13N+9,  13N+11,  13N+13 

17 

8 

17N+1 ,  17N+2,  17N+4,  17N+10,  17N+12,17N+13,17N+14,  17N+17 

19 

9 

19N+2.5 ,6 ,11 ,13, 15, 16, 17, 18 

23 

11 

23N+1, 4, 5, 8, 10, 15, 17, 18, 19, 21, 23 

29 

14 

29N+1, 2, 3, 4, 6, 10, 13, 17, 19, 20, 21  , 22, 25, 27 

31 

15 

31N+3, 5, 7, 8, 9, 14, 17, 21, 22, 24, 25, 27, 28, 29, 30 

37 

18 

37N+3, 5, 6, 7, 8, 12, 17, 21, 22, 23, 24, 26, 29, 30, 32, 34, 36, 37 

P 

(p-n/2 

PN+Mj  ,^2 ,M^ . . »M(p_i)/2 

♦Subscripts  such  as  1,2,3  etc.  in  3No+3,  SN^+l, 
elswhere  were  omitted  for  the  sake  of  clarity 


as  shown  in  Table  3  and 


congruent  to  0(mod  3).  If  just  the  divisor  3  were  to  apply  to  the  column  it 
still  could  only  have  a  maximum  Dc  of  66  2/3%  (65%  for  x=l  to  40).  The  latter 

set  applies  to  the  2nc*,  5^,  8^, .  columns;  here  two  of  every  three 

integers  of  a  column  are  divisible  by  3.  The  highest  density  possible  is  only 
33  1/3%;  because  of  the  presence  of  other  divisors  for  both  sets,  Dc  will 
usually  be  much  lower.  The  description  above  was  the  basis  of  the  criteria 
formulated  in  reference  1  which  indicated  which  equations  would  tend  to  be 
prime-rich.  The  criteria  were  that  c  be  prime  and  that  Xi,  the 
first  value  of  x  used  in  any  column  of  the  diagonal  array  of  Table  1, 
reference  1 ,  be  e  0  (mod  3).  The  latter  criterion  still  holds;  it  is  just 

another  way  of  stating  that  columns  3N+3,  N=0,l,2 .  with  no  integers 

=  0  (mod  3)  are  more  apt  to  be  prime-rich  than  columns  described  by  the  sets 
3N+1  and  3N+2.  The  second  criterion  still  applies  approximately  when  c  is 
relatively  small  but  as  c  becomes  appreciably  larger  it  is  composite  for  an 
increasing  proportion  of  prime-rich  equations.  Therefore  if  a  search  for 
prime-rich  equations  is  limited  to  values  of  D  >  67.5%  only  one  of  every  three 
columns  (columns  and  equations  are  used  interchangeably  here)  need  be 
exami ned. 

To  carry  the  above  analysis  one  step  further  one  can  consider 
coincidences  of  empty  columns  for  two  or  more  PCA's.  For  example  there  is  one 
set  of  empty  columns  for  P=3,  namely  3N+3  and  two  for  P  =  5,  5N+1  and  5N+4.  It 
is  easy  to  see  that  there  are  empty  columns  common  to  the  pair  3N+3  and  5N+1 

and  also  to  the  pair  3N+3  and  5N+4.  For  the  two  pairs,  for  N=0,l,2, . 

3N+3  =  3 ,6  ,9 ,1  2, 1  5 , 18  ,21 .  and  5N+1  =1  ,6 ,1 1  ,1 6  ,21 . .  The  equation  for 

columns  empty  of  both  the  divisors  3  and  5  is  given  by  ^  = 1 5 N + 6  and 
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TABLE  3.  COINCIDENCES  OF  EMPTY  COLUMNS  FOR  DIVISORS  3  AND  5;  AND  3,5,7;  AND 
3,5,7,11  FOR  DERIVED  PRIMITIVE  CELL  ARRAYS 


No.  of 
K  family 


Equations  for 
"empty"  columns  for 


3Nx+3 


5N2+1 


5Ng+4 


7N3+2 


7N3+6 


7N3+7 


7N3+2 


7N3+6 


7N3+7 


iin4+i 

"  4 

"  8 
"  9 

"  10 


iin4+i 


iin4+i 


"  4 

"  8 
"  9 

"  10 
11N4+1 
"  4 

"  8 
"  9 

"  10 


iin4+i 


"  4 

"  8 
"  9 

"  10 
iin4+i 

"  4 

"  8 
9 

10 


Families  of  equations 
of  form  x  -x+c:  c= 


2310N+2201 
"  "  1151 


"  "  521 

"  "  941 

"  "  1361 


"  "  221 
"  "  1481 

"  "  851 

“  "  1271 

"  "  1691 

"  "  881 
"  "  2141 

"  "  1511 

"  "  1931 

"  "  41 

"  "  1277 

"  "  227 

"  "  1907 

"  "  17 

"  "  437 

"  "  1607 

"  "  557 

"  "  2237 


m 


u‘i 

•5S!I 


8 


ss 


m 


A-'.’.SV 


v!W 
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similarly  for  the  second  pair  we  have  ^=15N+9,  N=0,l,2,...  .  As  c=2N3  ^-1 


have  the  equations  x  -x+30N+l<2»  N=0,l,2,..  where  I<2=11,  17.  Thus  there  two 


families  of  equation  for  coincidences  of  columns  for  which  there  are  no 


integers  congruent  to  0(mod  3,5).  For  convenience  this  is  called  a  (3,5) 


net.  Similarly  we  can  construct  the  net  (3,5,7)  for  coincidences  of  empty 


columns  for  divisors  3,5  and  7.  It  is  represented  by  I=x  -x+210N+Kg, 


N=0,l,2, .  This  represents  six  families  of  equations  where 


Kg=l 1 ,17,41 ,101 ,137,  and  167.  The  next  net  is  that  for  (3,5,7,11)  and  its 

O 

equation  is  I=x  -x+231 ON+K^g.  K^q  represents  30  families  of  equations.  These 


are  listed  in  Table  3.  This  can  be  continued  to  more  restrictive  nets  such  as 


the  next  in  succession,  (3,5,7,11,13)  whose  equation  is  x  -X+30030N+K, on.  But 


to  work  with  this  equation  and  succeeding  equations  one  must  evaluate  unwieldy 


numbers  of  constants.  For  example  for  K1Rn  there  are  180.  One  member  of  this 


O  c  4- 

family  is  Euler's  equation  I=x  -x+41  (this  is  the  215  column  of  the  original 

1  O 

array  of  Table  1,  ref.  1).  The  next  member  is  I=x  -x+30071  (N=l  in  the 


(3,.. ,13)  equation  just  above);  and  it  would  correspond  to  the  15036tn  column 


of  the  same  array.  Because  of  these  considerations  attention  was  focused 


primarily  on  the  family  I=x  -x+2310N+Kon  (and  to  a  lesser  extent,  on 


x2-x+210N+K6  and  x2-x+30N+K6 


CALCULATIONS 


Initially  was  held  to  the  value  40  and  N  was  varied  for  each  of  the 
30  values  of  K.  For  N=0  to  29  this  generated  900  equations  and  Dc  was 
calculated  for  each.  The  average  value  of  Dc  for  all  900  was  47.531;  for  1800 


v.v.v 


I'V.-.V 


(N=0  to  59)  it  was  about  41.72%.  Ultimately  the  range  for  N  was  extended  from 
0  to  650.  In  this  extended  range  the  largest  c  value  encountered  was  about 
1,500,000.  A  similar  evaluation  was  made  for  the  (3,5)  and  (3,5,7)  nets.  But 
the  largest  value  of  c  reached  for  b  h  was  only  about  84,000;  for  the  first 
net  N  was  =  3000  and  for  the  second  it  was  =  400. 

Most  of  these  calculations  were  accomplished  with  a  Basic  program  using 
an  HP9826  computer.  The  program  was  set  up  to  determine  the  values  x=l  to  40 
for  each  equation  and  included  a  Fast  Prime  Search  program  as  a  subroutine. 

The  latter  program  was  developed  by  S.  Varsano  and  can  determine  the 
primality  of  any  integer  less  than  100,000,000  with  double  precision.  However 
because  the  particular  computer  used  was  not  furnished  with  adaptation  for 
double  precision  the  values  of  Dc  (see  Table  1)  could  only  be  determined  for 
xN  values  up  to  about  xN=8,000.  The  results  for  x^=10,000  were  obtained  by  the 
use  of  the  Laboratory's  720  Cyber  computer  and  all  the  data  shown  in  Table  1 
were  confirmed  by  it.  After  calculating  values  of  I  and  indicating  whether 
prime  or  composite  for  several  hundred  equations  the  program  was  modified  to 
determine  the  percentage  of  I's  per  column  which  were  prime  rather  than  print 
out  each  value  and  label  it  prime  or  composite.  As  the  range  of  N  was 
increased  only  those  percentages  equal  to  or  greater  than  >  67.5  were  printed 
out.  Most  of  the  equations  having  high  Dc  values  were  detected  by  all  three 
"nets"  to  the  extent  that  equivalent  ranges  were  covered.  However  there  were 
exceptions.  For  example  the  (3,5)  net  picked  up  x  -x+107  (D^gy=75%)  which  the 
other  two  did  not  (and  could  not  by  definition  for  although  the  column 
represented  by  this  equation  does  contain  integers  divisible  by  7  and  17  it 
contains  no  integers  which  are  congruent  to  0(mod  3,5,11,13,19,23,37,41  and 
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others) ) . 

Restricting  attention  to  one  family  of  equations  alone  will  not  pick  up 
all  prime-rich  equations  but  will  find  most  of  them.  The  most  complete 
selection  is  assured  by  concentrating  on  the  (3,5)  net.  However  there  is  a 
tradeoff.  The  (3,5)  net  takes  considerably  more  time  (many  more  equations  are 
processed)  than  the  (3,5,7)  or  succeeding  nets. 

RESULTS 


Eighty  seven  equations  were  found  for  which  was  Dc  >  67.5%.  Of  course 
Euler's  equation  was  one.  The  next  highest  Dc  value  found  was  that  for 
I  =  x2-x+21  377  ;  D2 1 3 7 7  =  •  Four  other  equations  shared  the  next  highest  values 

of  Dc=77.5%;  sixteen  had  75%  etc.  Some  percentage  values  for  the  case  of  c=41 
are  well  known;  for  x=l  to  x^=100  0^  =  86%  and  for  x^=2398  D^^  =  50%  exactly.  In 
order  to  compare  values  of  Dc  with  that  for  Euler's  equation  x^  was  extended 
to  100,200, . ,10,000.  As  was  pointed  out  fifteen  of  the  equations  of 

O 

Table  1  were  more  prime-rich  than  x  -x+41  and  five  were,  to  a  close 
approximation,  almost  as  prime-rich  when  x^=10,000.  Published  results  (Ulam 

o  7 

et  al  )  state  that  D41=47.5%  for  I  =  10  ,  *N=  3162.  The  results  as  shown  in 
Table  1,  row  1  have  extended  this  to  0^  =  41.49%  for  I  =10^, x^=10 ,000. 

Figure  1  is  a  semilog  plot  of  the  density  D  ,  c=41,  21377,  etc.  as  a 

function  of  x^  (x^=  40,100, . 10,000).  The  dependent  variable  Dc  is  plotted 

as  the  abcissa.  Only  four  of  the  twenty  values  of  c  which  appear  in  Table  1 
are  graphed  in  Figure  1.  As  the  equations  I=x  -x+c  are  Diophantine  in  nature 
(the  I's  are  integers)  the  initial  parts  of  the  curves  where  xN=40 ,100, .. .are 
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FIGURE  1.  SEMILOG  PLOTS  OF  De  VS  XN  FOR  4  VALUES  OF  C  FROM  TABLE  1.  THE 

FIFTH.  SHORTER,  CURVE  IS  AN  AVERAGE  OF  ALL  20  VALUES  OF  THE  TABLE. 


10 
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only  very  rough  approximations  of  the  discontinuous  way  in  which  Dc  varies 
with  xN.  However  if  values  considerably  larger  than  xn=10,000  were  to  be 
examined  the  curves  although  still  discontinuous  would  become  much  better 
approximati ons  to  smoothly  varying  curves  and  would  probably  have  slopes 
approaching  each  other  in  value. 

ADDITIONAL  OBSERVATIONS 

2 

1.  C =4 1  occurs  in  every  family  of  equations  such  as  x  -x+6N+5  (N=6), 
x^-x+210N+41  (N=0)  and  on  up  to  x2-x+7 ,420 ,738 , 1 34 ,81 ON+41  (N=0)  the  last  of 

which  is  for  the  (3,5,7, . ,37)  net.  <2  is  the  constant  for  this  equation 

which  means  that  41  is  just  one  of  the  538,876,800  values  of  z! 

2.  For  the  last  named  family  of  equations  in  1.  above  the  second  member  is 

2 

obtained  by  setting  N=1  i.e.  I  =  x  -x  + 7  ,420 , 7  38 , 1  34 ,851 .  It  is  quite  possible 
that  for  x=l  to  40  no  integer  of  this  set  is  prime.  This  indicates  the 
impractical ity  of  a  direct  search  for  the  next  equation  in  the  above  family 
for  which  no  integer  in  the  evaluated  equation,  x  -x+41,  is  congruent  to  0(mod 

3,5,7, . ,37).  One  is  more  apt  to  find  such  an  equation  incidently  by 

using  nets  such  as  (3,51,  (3,5,7)  or  (3,5,7,11).  x  -x+107  was  found  in  this 
manner  (see  page  ).  Similarly  x  -x+21377  was  found  to  be  the  10Ln  member  of 
the  family  x^-x+231 UN+587  (N=9):  it  is  also  the  first  member  of 

x^-x+200,660,490,120N  +  Kz  which  is  the  equation  for  the  (3,5,7, . ,31)  net 

(that  is  for  N=0,  K=21377  and  z  represents  29,937,600  constants). 
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3.  Two  rules  for  calculation  were  found  with  regard  to  primitive  cells 
(by  primitive  cell  is  meant  the  area  of  the  primitive  cell  itself  comprised  of 
the  P  x  P  spaces  and  not  the  whole  primitive  cell  array).  These  rules,  though 
not  proved,  work.  The  first  is  that  the  number  of  empty  columns  in  a  cell  is 
equal  to  (P-l)/2.  The  second  is  that  the  number  of  coincidences  of  empty 

columns  is  given  by  (P^-l )/2  x . x  ( ) /2  x  (P^-l )/2.  For  example  for 

PN=11  Nec= ( 1 1 - 1 ) /2  x  (7-1 ) /2  x  (5-1 ) /2  x  (3-l)/2=30.  This  gives  the  number  of 

O 

values  of  the  constant  K  required  to  be  evaluated  in  x  -x+2310N+I<3q  (see  the 
last  column  of  Table  3). 

P=P 

p 

4.  The  general  equation  for  any  net  is  x  -x+N  n  P+K^;  K(p)  is  the 

P=2 

number  hut  not  the  values  of  the  constants  to  be  evaluated  and  is  given  by 

( Pp-1 )  /2x . . ,x  (P3-I  )/2  x  (P2-I  )/2.  Special  cases  are:  if  P=2  then  -1 

and  if  P=3  thep  k(1)  =  ^  (expressions  P^  in  3.  and  Pp  here  are  i nterchangabl  e) . 
For  P=  1 3  I  =  x^-x+N(  1  3x1 1  x7x5x3x2)+K^ =x^-x+30030N+K^g0  where  is  equal  to 

(13-l)/2  x  (ll-D/2  x  ( 7-1  )/2  x  (5-1 ) / 2  x  (3-1  )/2  =  180  i.e.  K ( p )  =K ( i 3 )  =Ki80 

5.  It  is  worth  noting  that  for  all  87  equations  having  values  of 

Dc  >  67.5?,  20  of  which  are  included  in  Table  1,  the  last  digit  for  every 
constant  c  is  either  a  1  or  a  7.  This  is  a  consequence  of  how  the  nets  are 
used  to  find  prime-rich  equations.  For  the  (3)  net  i.e.  for  p=3  in  4  above, 

p 

I =x  -x+6N+5 ,  N=0,l,2....  can  have  any  odd  integer  as  a  unit  digit.  However  in 
order  to  more  readily  determine  which  column  (equation)  would  be  prime-rich  it 
was  necessary  to  go  to  the  (3,5)  and  succeeding  nets.  For  the  (3,5)  net 


1  2 
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I=x  -x+30N+ll,17.  For  the  two  families  represented  here,  where  N=0,l,2 . . 

every  c  value  will  be  an  integer  whose  units  digit  is  1  or  7.  As  all  succeeding 
nets  involve  multiples  of  the  repeat  interval  for  coincident  empty  columns  for 
the  (3,5)  net  any  c  value  for  all  nets  beyond  the  (3,5)  net  will  also  have  a  1  or 
7  for  the  units  digit. 


CONCLUSIONS 

It  has  been  shown  that  there  are  only  six  equations  of  the  type 
2  2 

I=x  -x+c,  x=l,2,3,....(x  +x+c,  x=0,l,2 .  is  an  equivalent  form)  for  which 

all  values  of  x  from  1  to  c-1  yield  prime  values  of  1.^  They  are 
c=2,3,5,ll ,17  and  41.  c  =  l  could  be  a  seventh  but  the  solution  I  =  f ( x )  =  1  is  not 
prime.  Of  the  above  Euler's  has  been  the  most  prime-rich  equation  known. 
However  the  above  constraint  does  not  rule  out  the  possibility  that  other 
equations  of  this  form  can  be  just  as  rich  or  more  so.  In  fact  it  has  been 
shown  above  (Table  1)  that  there  are  Such  equations  and  that  there  is  a  viable 
method  of  finding  them.  Work  is  being  completed  on  an  analogous  set  of  at 

•p 

least  16  equations  of  the  related  type  I=x  +x-r  e.g.  for  r=501,229  and 
x^=10,000  Dp  =49.25%.  This  equation  is  a  "row"  equation  of  the  special  array 
described  in  reference  1  in  contrast  to  the  "column"  equations  discussed  in 
this  report. 

It  should  be  feasible  to  examine  c  values  greater  than  1.5  x  10^. 

Between  c=858,707  of  Table  1  and  c=  1.5  x  10^  the  only  c  value  found  for  which 
Dc  was  equal  to  67.5%  was  c=l,275,707;  Dc  for  x^  =10,000  was  about  37.5%.  It 
is  anticipated  that  there  will  be  few  if  any  Dc's  J  67.5%  between  1.5  x  106, 
the  highest  c  value  examined,  and  3.0  x  10^. 

It  should  also  be  worthwhile  to  extend  the  range  of  xN  beyond  10,000 
for  all  values  of  c  in  Table  3.  Using  only  the  five  viable  points  available 


13 


for  for  x^=l,000  to  10,000  the  following  empirical  equation  was  found: 

j  =  141.9(  xN'°*  .  The  value  of  for  x^=ll,000  has  been  found  to  be 

40.96..%  ;  that  is  4506  out  of  11,000  I's  (integers)  are  prime. ^  Substitution 
of  xN=  11,000  in  the  above  equation  yields  D4j=40.78  (or  4486  primes  out  of 
11,000  I's),  which  differs  from  the  true  value  by  less  than  0.5%. 

Substitution  of  x^=  100,000  in  the  same  equation  yields  =  30%.  However,  how 
valid  this  figure  may  be  depends  on  obtaining  more  points  in  the  relation 
Dc=f(xN)  and  finding  a  better  fit  to  the  plotted  curve. 

Two  other  items  of  possible  interest  related  to  this  report  which  may  be 
worth  looking  into  are  to  examine  those  columns  and  rows  (equations)  which  are 
very  prime-poor  and  to  follow-up  a  number  of  equations  which  were  found  to  be 
prime-rich  but  fell  just  short  of  the  arbitrarily  set  Dc  criterion  of  being 
>  67.5*-.  One  such  equation  was  for  c  =  l,384  ,907.  The  Dc  values  for  x^  =  40, 
2,398  and  8,000  were  65,  42.66..  and  38.725%  respectively  and  the  projected 
density  for  x^=10,000  is  about  36%. 

*  The  Cyber  720  gave  40.97....%;  for  x^=12,000  it  gave  40.35%. 
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